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AlisTRACT 


It is necessary to have an understanding of the charge-exchange plasma 
environment around a spacecraft that uses mercury ion thrusters for propulsion 
so that the interactions between this environment and the spacecraft can be 
determined. A model is reviewed which describes the propagation of the mercury 
charge-exchange plasma and extended to describe the flow of the molybdenum com- 
ponent of the charge-exchange plasma. The uncertainties in the models for vari- 
ous conditions are discussed throughout this report. Such topics as current 
drain to the solar array, charge-exchange plasma material deposition, and the 
effects of space plasma on the charge-exchange plasma propagation are addressed. 
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A. INTRODUCTION 


Ion thrusters will be used in the near future as a primary propulsion device 
for interplanetary spacecraft and for attitude control and orbit transfer of 
spacecraft near earth. Because of the plasma environment produced by the thruster 
it is important to understand how this propulsion device will affect various 
spacecraft subsystems. An interplanetary spacecraft will likely be most affected 
by the ion thruster environment because of the use of multiple thrusters and 
thrusting periods of years. 

This paper will describe the charge-exchange plasma environment produced by 
the thrust subsystem and how this environment will interact with the spacecraft. 
The effect that the space plasma will have on the charge-exchange plasma will 
also be addressed. 

In order to describe the thruster produced environment, the various ion 

species emitted by a thruster are described. To distinguish the various ion 

species emitted by an ion thruster the notation of group I, group II, group III, 

and group IV ions was introduced.^ The group I ions are those of the primary 

beam which travel dovmstreara within a cone that has a small divergence. They are 

accelerated ty the total potential difference between the bombardment discharge 

and the neutralized thrust beam. Group II ions are energetic ions that leave a 

thruster at angles up to nearly 90 degrees to the beam axis.^ They are formed by 

charge-exchange between the screen and accelerator grids of the thruster optics. 

The group III ions are charge-exchange ions that are formed just downstream of 

the accelerator grid. These ions are draTO back into the acceleration grid by 

its negative potential. This produces a sputtered efflux from the accelerator 

optics that is somewhat more peaked (in the downstream direction) than a cosine 
2-5 

distribution. The neutrals which leave the thruster do so in approximately a 

g 

cosine distribution. All of the preceding effluxes are essentially limited to 
llne-of-slght trajectories from the ion optics of the thruster. 

Two other efflux categories, the group IV ions, are not limited to line-of- 
sight trajectories, and therefore can pose special spacecraft contamination 
problems. These effluxes are the low-energy mercury and molybdenum charge- 
exchange ions that are generated downstream of a thruster. Because these 
effluxes are low energy ions, the trajectories can depart substantially from 
straight lines due to only moderate electric fields. The natural tendency is 
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for less dense regions of the charge-exchange plasma to be negative relative to 
the dense regions, so that ion trajectories are bent from the more dense to the 
less dense regions. l.s a result, the entire volume around a spacecraft becomes 
filled with charge-exchange plasma, with the density of this plasma depending 
on the ease of access from the region of generation downstream of a thruster. 

The charge-exchange plasma is mccroscoplcally neutral so that everywhere electron 
density equals ion density. 

The propagation of a charge-exchange plasma from a thruster has been the 

S'" xo 

subject of study due to the possible Interactions with spacecraft. Further 
study of this propagation Is expected. This report is Intended as a summary of 
the present understanding of this propagation, interaction effects thereof, and 
possible verification experiments. The Jimitations of this present understanding 
will also be given, where possible. The calculations presented in this re 3 >ort 
will be limited to mercury propellant and molybdenum grid material, although the 
process with other materials would be expected to differ mostly due to different 
cross sections and yield coefficients. All equations presented herein are in 
SI (mks) units unless otherwise stated. 

B. MERCURY CHARGE-EXCHANGE MODEL 

The llg charge-exchange model presented here is essentially the same as 
one presented in a previous publication. The assumed point of origin for 
the Hg charge-exchange plasma is Indicated in Fig. 1. For a single thruster, 
this point is one thruster radius, rj;, downstream of the ion optics and on the 
thruster centerline, as indicated in Fig. 1(a). This point was selected because 
it gave the best correlation of experimental data from several tests. 

Ft>r a circular array of thrusters, the effective point of origin would be 
expected to be one radius of the thruster array, r^' (see Fig, 1(b)), doVmstream 
of the center of the array. For a noncircular array, the effective origin would 
be expected to be downstream of the center of the array, between half the maxi- 
mum array dimension and half the minimum array dimension. A noncircular array 
would also be expected to cause some departures from axial symmetry in the dis- 
tribution of Hg charge-exchange plasma. But, because Hg charge-exchange measure- 
ments have not been made with thruster arrays, the exact extent of these 
departures from axial symmetry is not known. 
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Figure 1. Parameters for Calculation of Charge-Exchange 
Plasma Generated by Thruster. 
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To calculate the electron density and saturation electron current density, '■ 

select parameter values from Tables I and H for the angle 9 (see Fig. 1) of 

. i ; 

Interest and substitute in the equations bfvlow: ' 

"e " - Huy^b^^nuHTable I Param.], m"^ (1) , 

^e " n,i>/rb^l\jlTable II Param.], Am~^ (2) 


Table I. Electron Density Parameter as a Function of Angle 
from a Beam Direction 


e. Deg. 

Parameter 

0, Deg. 

Parameter 

90 

12 

2.5 X 10 

140 

11 

3.2 X 10 

100 

I"? 

1.8 X 10 

150 

1.9 X 10^^ 

110 

12 

1.2 X 10 

160 

1.1 X 10^^ 

120 

12 

8.2 X 10 

170 

6.7 X 10^® 

130 

11 

5.1 X 10“ 

180 

3.8 X 10^° 

Table II. 

Electron Current Density 
From Beam Direction 

Parameter as 

a Function of Angle 

0, Deg. 

Parameter 

0, Deg. 

Parameter 

90 

0.12 

140 

0.017 

100 

0.091 

150 

0.011 

110 

0.064 

160 

0.0059 

120 

0.044 

170 

0.0034 

130 

0.028 

180 

0.0022 


The beam current, is in amperes (A). It is the total current such that 
if a beam is produced by six thrusters each operating at 2A, J^, will be 12A. The 
propellant utilization, is a dimensionless quantity which can be expected to 
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Jk- “0.9. The beam radius, r^, and the radius R are defined by Fig. 1. Both are 
measured in meters. 

Tables I and II represent the mean values of experimental data. The 
upstream hemisphere (90-180 deg. from beam direction) i'pt the region of most 
interest. In this hemisphere the data tanged upwards and downwards from the 
values used for Tables I and II by about a factor of two. 

For electron density, the electron temperature i.^ also a consideration. 

The production rate of charge-exchange ions is independent of electron tempera- 
ture, but the velocity of departure for these ions varies as the square-root of 
electron temperature. To the first approximation, then, ion (and electron) 
density will vary inversely as the square-root :al~etroa temperature. From 
available data, the electron temperature i:i the charge-exchange plasma around the 
thruster is about half of the electron temperature in the ion beam. The data 
us6m for Tables I and II had electron temperatures in the charge-exchange plasma 
of 2.5 and 3.5 eV. These values are believed representative of 5-15 cm Hg Ion 

thrusters. The electron temperature within a 30-cm Hg thruster beam was found 
11 

to be 0.35 eV. Including the factor of two lower electron temperature in the 
charge-exchange plasma, a 30-cm thruster would be expected to have about a 
factor of four higher electron density than given by Eq. (1) with the parameter 
values listed in Table I. 

The scaling of electron density inversely as the square-rcoc of electron 
temperature, as described above, is actually of limited accuracy. Another con- 
sideration is the directed neutral velocity before charge exchange. Because the 
neutral temperature varies only slightly between different sls'.es of thrusters, 
this initial velocity is much more significant where electron temperatures are 
small. We do not know the magnitude of this effect, but we do know its direc- 
tion. The initial downstream neutral velocity of charge-exchange ions will bias 
their distribution towards the downstream direction, thus reducing their density 
in the upstream hemisphere. Because of the smaller plasma electron energy in 
large thrusters, this effect will be most noticeable in such thrusters. The 
additional effect of neutral velocity, beyond the density varying inversely as 
the square-root of electron temperature, should make the actual electron density 
lower than the calculated values for 30-cm thrusters in the upstream hemisphere. 
At present, the magnitude of the reduction is not known. 


Electron CoraperaCure Is not as Important for t:he saturation electron current 
density. Except for the bias duo to initial neutral velocity, the effects of 
electron temperature cancel for electron current density. 

The ion arrival race is closely related to Che saturation electron current. 

The ion velocity can be closely approximated by the Bohm, or ion acoustic, veloc- 

8 Q 

ity based on electron temper iture in the ion beam. ’ The electron velocity is, 
of course, a Maxwellian distribution corresponding to the temperature in the charge- 
exchange plasma. For a surface normal to the local ion velocity, the ratio of 
ion current to electron saturation current is 


- 2 = 5.9 x 10“^. (3) 

It sliould be kept in mind that tht ion velocity is substantially directed, 
while that of the electrons I'a isotropic. With tliis limitation in mind, Eq. (3) 
can be used to obtain ion current detusities from Eq, (2). For convenience the 
following extension to Eq. (3) gives the time (t) in seconds for a monolayer 
to form. This is assuming that the surface is cold enough to condense all the 
mercury arriving there. Also, the tlieory does not allow for geometrical con- 
siderations such as shielding. Shielding may act to reduce Che arrival rate 
but the charge-exchange plasma will be able to flow around obstructions, 

t - 3.3 X 10^/j^ (3a) 

1. Sample Calculation for Spacecraft 

A sample calculation, using the charge-exchange plasma model described, is 

performed to illustrate the variation in plasma and current densities around a 

typical spacecraft. The spacecraft has an assumed configuration which utilizes 
'12 

eight 30-cm thrusters arranged in two parallel rows of four each. Each row 
of four has center-to-center spacings of 76 era, while the centet-to-center spacing 
between the two rows of four is 51 cm. The two solar array panels extend 
radially outwards from the axis of Che spacecraft, with the inside end of the 
panels starting 254 cm from the spacecraft centerline and 432 cm upstream of 
the thruster exhaust plane. 


6 


The calculation procedure used is that described above. Seven of. the eight 
thrusters are assumed to he operating at a beam current of 2A each with a 
propellant utilization of 0.9. The other thruster is not operating. Although 
the nominal diameter is 30 cm for these thrusters, the actual beam diameter is 
closer to 28 cm. 

The charge-exchange plasma is assumed to have an electron temperature 
equal to half the measured ion-beam electron temperature (from single thruster 
tests) of 0.35 eV. This lower electron temperature, compared to the data used 
for Tables I and II, should result in an electron density increase of about 20*. 
As mentioned above, there is no similar correction for the current density 
parameter. Another assumption concerns the point of origin for the charge- 
exchange plasma. 

For the non-circular thrust-^r array of this spacecraft, the totjjl array 
area is equivalent to a circle wi;h a radius of 0.885 m. The effective origin of 
the charge-exchange plasma is assumed to be 0.885 m downstream of the center of 
the thruster array. Some nonsymmetry of the charge-exchange plasma about the 
thrust axis is expected to result from the nonsymmetrlcal distribution of 
thrusters, but the charge-exchange plasma will tend to redistribute itself 
so os to reduce this effect. The maximum increase in electron density or 
electron current is therefore estimated at less than a factor of two for a 
distance of one meter or more from the thruster array. (Th^s factor increases 
the over.Tll uncertainty to about a factor of four.) 

The plasma properties at the centerline of the solar array were calculated 

as a function of distance from the inboard end of the solar array and plotted in 

Fig. 2. The combined effects of radial distance, R, and angle, G, give maximum 

electron density and current density near, but not at, the inboard end of each 

solar array. For the assumed dimensions of 3.8 ni x 44,3 m for each array, the 

integrated total of the current density shown in Fig. 3 is 1.7 A fbr both arrays. 

This current includes no increase In effective array size due to sheath thidfe- 
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ness. The electron density ranges from about 1 to 12 x 10 m , viilch, at 
0.175 eV, corresponds to a range from 3 mm to 1 cm in Debye distance. For a 
solar array that is 400 V positive of the surrounding plasma, for example, the 
total she.ith thickness would range from about 20 co 60 cm. The effect of sheath 
thickness can clearly be significant, but will not change the qualitative effects 
shown in Fig. 3. 





Tlie locations on tlu' spacocraft correspond to wide ranges of both radial 

distance, R, and angle from beam direction, 0, Two dlnionslonal plots of electron 

density and saturation electron current density are slnn^i In Figs. '♦ and 5, Holli 

plan and aide views of the assumed .spacecraft '.ue shown In each of those figures. 

The as.sumed point of origin for tlie charge-ewclumge plasma is Indicated In Figs. 

and 5 and is, a.s was mentioned previously, OtBRS m dewnstroam of the center of 

the thru.ster-array exhaust plane. As shown in Fig. /», the expected electron 

'i 7 -■} 

density can range from about in to 10 cm . In Fig. '5, the saturation electron 

9 

currant density is shown to range from 1 to 100 iiiA/nr. From Fig. 5 and Rq, (3) 

2 

the corresponding arrival rate of llg ions would range from 6 to 600 pA/m . Fig- 
ure 6 shows the arrival rale of Hg In terms of monolayer formation. One 

0 

monolayer of mercury is 'v2.9A. 


As mentioned, seven operating thrusters were assumed for the calculations 
whose results are showi in Figs. 2 through 6. For fewer thrusters or less than 
rated operating conditions, both the eleetron densities and the electron and ion 
current densities would have to be re<lucod. For fewer operating tlirustors, 
these parameters should be reduced in nroportlon to the nvimber of thrusters. 

For less than rated eondJtions, the plasma parameters would he proportional to 
For the usual throttled operation m-ar the dlscliarge- 
perfornianeo "knee", will bt> nearly proportional to 

The orientation of a surface relative to tiie propagation of the charge- 
exchange plasma will Influence the density and arrival rates observed at tlie 
surface. Some reduet isn of these parameters would be expeeteii for more pro- 
tected surface or ientatlons and surfaces that are locally shiolded by parts of 
the spacecraft. Ileeause of the proven capability of a charge-exchange plasma 
to flow around corners, the amount of this reduction cannot be estimated at 
die present time. 

A variety of devices lias been suggested and tried for diverting or collect- 
ing the charge-exchange plasma. Shields near tiie thruster cun dis.p lace or divert 

tlie charge oxchange plasma, but do not appear to significantly reduce the number 

' 9 

of cliarge-exclinnge ions that escape to the surrounding volume, presumthly bet auso 
electric fields do not extend far into the plasma. Fositive surfaces will repel 
chtirge-excliange ions, but the large electron currents also collpcted will usually 
make tills approach impractical, Insulated surfaces will draw equal numbers of 
electrons and ions, which means that all of the arriving Ions will be collected 
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and most of the electrons will be reflected. Negative surfaceH will reflect 

electrons and collect only ions, but this effect appears to be of little value 
in significantly reducing the upstream flow of the charge-exchange plasma, 

2. Discussion of Sample Calculation 

The amount of electron current collected by a positive and exposed solar- 
cell array is one of the significant results of the calculation presented above. 
A high voltage array would be more seriously affected by a fixed level of 
electron current collection, so the results are also voltage dependent. For a 
400 V array, for example, the 1.7 A electron current would correspond to a maxi- 
mum loss of '^680 watts. If there is a significant variation in voltage along 
the length of the array, this loss can be reduced by placing the high voltage 
end of the array at the outboard end, where the electron current density is 
smallest. This approximation assumes there are exposed conducting paths for 
the electrons to be drawn from the charge-exchange plasma. 

2 

The arrival rate of electrons is a maximum of 14 niA/m on the solar array 
(see Fig. 3). By Eq. (3a), this arrival rate corresponds to a monolayer buildup 

O 

of mercury at Chat position In about 6 hoiirs. A 1000 A layer with near zero 
transmittance will be deposited in about 83 days. Significant effects are 

O 

likely to occur with layers less than 1000 A. From experience with mercury 
vapor, we know that no accumulation occurs when the surface involved is above 
the evaporation temperature for that rate. There are still two reasons to be 
concerned though. If the surface should fall below this temperature, an 
accumulation will occur that will require a considerable excess of temperature 
to evaporate again. This could very likely be the case for the exterior of the 
sun-shaded sides of thermal blankets, passive radiators, and other spacecraft 
surfaces. Onc^ the mercury is accumulated, it may be difficult to evaporate 
away. This is because mercury evaporation can be orders of magnitude below 
equilibrium values when the surface is less than very cl.e^n and the temperature 
is below roughly 100 “C. 

A less clear aspect of possible Hg accumulation is that the arriving Hg is 
ionized and is probably arriving with several eV of energy. Because 1 eV corre- 
sponds to 11,600°K, we siiould not be too hasty in assuming that nc damage will 
occur. 
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3. Llnti cations of Model Accuracy 

Most of tlie limitations associated with the Hg charge-exchange model and 
calculation have been mentioned earlier. They will be summarized here for 
completeness . 

The model was obtained from tests of rather simple configurations, hence 
is not suited to estimating ths effects of nonsvmmetrlcal thruster arrays and 
complicated spacecraft shapes. Even the effects of surface orientation are not 
fully understood. 

The electron temperature it the charge-exciiange plasma is an important, but 
only partially evaluated pararaecsr. We know that smaller thrusters tend to have 
higher temperatures in this plasna, but we do not know if arrays of thrusters 

will depart from the values obtained with single thrusters. We also know 

* 

that initial neutral velocity v/lll be more important at small electron tempera- 
tures, resulting in a bias of chai'ge-exchange ion trajectories in the downstream 
direction. The magnitude of this bias, though, has not been evaluated. 

Finally, there Is an electron temperature effect that has not been mentioned. 
The space environment includes a plasma. The electrons of this space plasma 
will interact with the electrons of the charge-exchange plasma, possibly result- 
ing in a different electron temperature than obtained in ground tests. This 
effect will be discussed in the section, Electron Temperature in Space. 

The uncertainties that now exist in the model are summarized in Fig. 7. 

C. MOLYBDENUM CHARGE-EXCHANGE MODEL 

As mentioned in the introduction, a small fraction of the Mo atoms that 
are sputtered from the accelerator grid charge exchange with beam Hg ions down- 
stream of the thruster. These Mo charge-exchange ions share the ability of Hg 
charge-exchange ions to flow around obstacles. As will be shown, the rate 
of production of Mo ions is very low compared to that of Hg ions. It is there- 
fore difficult to make an experimental evaluation of Mo Ions. At this time, only 
a simple theoretical analysis can be presented, based on the Hg charge-exchange 
model presented previously. 
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The approach used la to compare the production rates of Hg and Mo charge- 
exchange ions. The Mo ions are assumed to be distributed in the. same manner as 
the Hg ions, but lower by a factor corresponding to the lower production rate. 
This approach Ignores the much higher initial velocities of the Mo neutrals. 
These higher velocities will result in a bias of trajectories in the downstream 
directiop. The arrival rates of Mo ions in the upstream hemisphere should, due 
to tills effect, be less than calculated. 

The production rate of Hg charge-exchange ions is 


N 

ce 


2Jb“(i - n j g 


(4) 


where q is the electronic charge, v^ is the average neutral velocity before 
charge exchange, a is the charge-exchange cross section, and ciio other param- 
eters are the sante as defined earlier. A similar production rate for Mo charge 
exchange ions is 


N 


ce,Mo 


2F J, h 
aba 

2 _ 


(5) 


where, in addition to previously defined variables, F is the ratio of accelerator 

<l 

impingement current to beam current, and Y is the sputter yield in atoms/itm. 

The ratio of Mo to Hg charge-exchange rates is obtained by dividing 
Eq. (3) by Eq. (4). With the use of Hg and Mo subscripts where appropriate, 
this ratio is 


N w 
ce,Mo 


K 


c«,Hg 


F Yn „ 

a u Mo o,Hg 


- 1u> “hb''o,Mo 


( 6 ) 


For the 30-cm thruster, the value of F^ is about 0,0025, “the value of Y for 
1000 eV Hg ions is approximately 1.0,^^ + Mo Hg + Mo^) is about 

1 X lO”"^^ + Hg Hg + Hg'^) is about 6 x 10 %,Mo 


is 


3 18 

abput 7 Jt 10 m/soc, '' ii about 230 m/sec (500"K) , and r) is again 0.9. 

o y (I g u 

With these substitutions, Eq, (6) becomes 


(7) 


The uncertainr; V in this calculation is due to the uncertainty of distribu- 
tion and yield of sputtered Mo for an operating thruster, the uncertainty of 
mean neutral velocity (the datum used here was for Kr on Mo), the uncert.iindy 
in bombarding ion energy, and the additional bias on Mo ion trajectories due 
to the initial high neutral energies. Ml of these additievojl uncertainties 
beyond those of charge-exchange Hg are estimated to raise the factor of uncer- 
tainty from four for Hg ions to ten for Mo ions. 

Equations (3) and (/), together with the electronic charge, c.in be used to 
translate the electron current density of Fig. i5) to Mo Ion arrival rate per 
unit area. The ratio is 


M,, /A 
Mo 


s . A X 10 


1,2 


Je* 


Mo 


-2 

m - sec 


( 8 ) 


This ratio was used to obtain Fig. 8. 
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A monolayer (2.5 A) of Mo atoms is about 1.6 x 10 /m* (as compared to about 
1.2 X 10^^/m^ (2.9 A) for Hg) . For the region of the solar array with the 
maximum of lA mA/m^ of electron current, then, a monolayer would be expected to 
take 72,000 hours. Even with a factor of 10 increase in rate due to the esti- 
mated inaccuracy of the calculation.^, it would take over a year for a monolayer 
to form. 

In situ measurements of transmittance have been made using films of the 
sputtered efflux from thrusters (mostly Mo). These 0,iiins were primarily from 
deposition of neutrals, rather than the charge-exchange Mo ions considered in 
this section. These results, however, should give an indication of the deposl- 
itlon effects to be expected. rhe total transmittance was not linear with time. 

It ranged from about 0.75 at 30 A thickness, to about 0,50 at 100 X, to about 

O 

0 at 1000 A. It appears, though, that assuming a linear variation Is conserva- 
tive below 30 % thickness. A monolayer of Mo should give about a 2,5 percent 
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reduction in transmittance. 




Arrival Rate, No'm -sec 


For the worst case (ten times the predicted rate) at the worst location on 
the solar array, over 1 year of operation with 7 thrusters would be required to 
Rive a 2.5 percent loss in transmittance. Because other locations should be 
less affected, it appears that No deposition on the solar array will not be a 
significant problem for the configuration studied, A much closer location to 
Che thru.ster, though, could cause problems. 

The major crnslderation for Mo deposition (in the upstream hemisphere) 

appears to bo sensitive instruments, optical or otherwise, that will be adversely 

affected by 10 to 100 percent of a monolayer. A monolayer on solar cell surfaces 
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with a normal emissivity (e) of 0.85 will produce e '^'0.7. Molybdenum deposi- 
tion could have an effect on thermal surfaces such as passive radiators used to 
cool ins tr ament component'.. 

D. CHARGE-EXCHANGE MEASUREMENTS TN GROUN"/ TESTS 
1. Mercury 

The Hg Ions -.a/kG up the bulk of the charge-exchange plasma. Measuremorits 
of this plasma density (with Langmuir probes) will therefcio give a close indica- 
tion of the density of Hg ions. As described in the experimental studies,^ 
the ions should be collected by negative boundaries, so reflection will not be 
a limitation on measurements, In fact, plasmas generally assume a potential 
more positive than surrounding surfaces, so that the proper negative wall bias 
is automatically attained. The only exception is for regions of very low density 
so that the total electron currents to boundaries associated with such regions 
is still quite small. Under these conditions the local plasma potential will be 
governed by wall potential. A moderate additional vvall bias (relative to the 
beam target) will avoid problems in the very low density regions. 

Some additional calculations or measurements arc desirable to make sure 
that most of the charge-exchange ions are generated immediately downstream of 
the thruster(s), rather than from the background neutral density all along the 
beam(s). To meet this latter requirement simply requires a low enough facility 
background pressure. 

The '"oasurement of Hg chargo-excliangc effects is thus readily accomplisiied 
in ground tests. 
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2. Molylidenum 

The measurement of Mo charge exchange ions poses much more difficult 
problems. One major problem Is the background level of sputtered Mo atoms, 

To show the magnitude of this problem, it is necessary to go through a simple 
analysis. 

The loss rate of sputtered Mo atoms from the accelerator grid is 

where all the parameters are defined the same as previously* This rate is 
probably an upper limit for No atoms in the doi^mstream direction in that not all 
accelerator grid Impingement is on the downstream side. 

Assuming a frozen Hg target in the ground facility to minimize contamination 
with other sputtered material, the arrival rate of Mo atoms at this target is 
due to an approximately cosine distribution from the accalerator grid. Sputter 
pits will develiop (midway between each group of three adjacent holes) on the 
downstream side of the accelerator grid. As this occurs, the distribution will 
depart from a cosine, but the departure will primarily be restricted to large 
angles from the beam direction. The rate of Mo arriving at the target shQald 
be nearly unchanged by these pits. The cosine distribution corresponds to a 
uniform arrival rate at the sphere shown in Fig. 9. The effective diameter of 
the beam at the target is d^., with a corresponding area of A^,. This diameter 
it? best determined experimentally, using the diameter over which the frozen Hg 
experiences net sputtering. Outside of this diameter the sputtered Mo will tend 
to be trapped within the condensing Hg vapor. From geometrical considerations, 
(Fig. 9) , the fraction of sputtered Mo atoms that will be resputtered from the 
target is approximately the solid angle Aj./(A^/2) divided by 
solid angle of Att. The resputtered fraction is thus 

F = A^/'iT(!^^. (10) 

rs t' t 



tal sphere 
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A similar fraction can be determined for the fraction of resputcered Mo atoms 
that arrive at a sample coupon. Assuming that this coupon is in the plane of 
the accelerator system and at a radius R from the thruster, with 

. ' (U) 

Multiplying Eqs. C9) through (11) to obtain the arrival rate of neutral atoms 
at che coupon area A^, 


N - F J.YA A ,/irV ^ • 
o,c ab ct ^t 


( 12 ) 


This rate can be compared to the arrival rate of charge-exchange Mo iops at the 
same location, using the preceding Mo model. From Eqs. (2), (3), and (6), 
togetlier with the Table II parameter for 90 de . , 


^ “ 7.1 X 10 ^ A F J “YOw V „ /qr,R^O,.v ,, 

ce,c can Mo o,Hg ^ b i.g o ,Mo 


(13) 


The neutral Mo to charge-exchange Mo arrival rate ratio is Eq. (12) divided by 
Eq. (13),. 


N /N = 140 A,.r, v a., v „ 

o,c ce,c Lb Hg o,Mo t b Mo o,Hg 


(14) 


Assume a 30-cm thruster, so Chat is 2A and is 0.14 ra. For the other values 

leg.), 0.5 m for R, and pr< 

With these substitutions. 


we will assume 3.5 m for 3.6 ra“ for A^. (±15 deg.), 0.5 m for R, and previously 

used values for q> o , c , v , „ , and y . 

^ Hg’ Mo o Hg’ ''o,Mo 


Eq. (14) becomes 


N /N = 21. 
o,c ce,c 


(15) 
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The Mo neutrals come from the target from tlie general downstream direction, 
while rhe Mo chnrge-exchange Ions are generally moving radially outwards, so 
some separation can be made by using the proper angle and slilcldlng for the 
coupon area, A^. With a "noise" level 21 times ulie "signal" level, thougli, it 
tpay be still difficult to separate out tlie Mo charge-exchangf tate by deposition, 
Also, although resputterod Mo auiin-v will have a "sticking" coefficient near 
unity, it will still be less tfian unity. A small percentage of Mo atoms reflected 
from the walls of the vacuum facility will add a nearly isotropic component th/it 
will bo hard to shield. 

Detection of Mo ions is a possible alternative, but: this will require 

4 

momentum analysis to discriminate from the nearly 10 times higlier background 
of Hg charge-exchange ions (Kq. (7)). Also the rate of Mo ions is small in 
absolute terms. For the location used above for comparison of neutral and 

9 

charged Mo, the current Mo ions throilglh '-1 cm" aperture (which is large for 
most momentum analysis) is only 10 A. 

In summary, then, experimental detection of Mo charge-e.xchnnge ions is 
difficult in a ground teat. 


E. ELECTRON TEMPERATURE IN SPACE 

The cliarge-exchange plasma in a ground facility fills, with varying 
degrees of density, the entire facility volume. The electrons in this plasma 
have long moan free patlis, so that wltli sufficient energy they can pass through 
any part of the available volume. The electrons in both the beam and the 
cliarge-exchange plasma travel much faster tlian the ions. Current neutrality 
v^ill therefore require, on the average, a number of electrons to be reflected 
at tiie boundaries for each one that escapes. 

In space, the charge-exchange plasma must evenctiaily blend Into the ambient 
space plasma, rather than be terminated abruptly at facility wails. A prelimi- 
nary model of the interactions or ion drive propulsion with this space plasma 
is presented in this section. 
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1. Space Plasma 


Two regions in space are considcrotl for this interaction with the ambient 
plasma. One is a low earth orbit, where an ion drive vehicle might be expected 
to start an orbit raising operation. The other is Interplanetary space at the 
radius of earth's orbit from the son. These two reg,ions are by no meann a 
complete inventory of ion drive environment. They are, however, surf'^cient to 
shew the major problems involved. 

The low earth orbit (MSO nautical miles) has low energy electrons (0.1 - 

6 **3 20 

0.2 eV) at a iilgh density ('^10“ cm ). For interplanetary space, the electron 

—3 2^ 22 

energy is much higher (10 - 20 eV) and the density much lower (5 - 20 cm ) . ' 


2. Distance of Interaction 


\ Tlio solar wind in interplanetary space travels radially outward fr'>m the 

Xa pi 

Sy\n aud the magnetic field lines are frozen into this plasma. ' The magnetic 

Re^yiolds number can be used to assess the interaction of this magnetic field 

with tite charge-exchange plasma. The magnetic Reynolds number is given by, 

R « h a V fi, , 
m 


where a is plasma conductivity, and v and Z are the typical velocity and lengtli, 

respectively, of the situation In que.stion. Taking the permeability, p = 

and for reasonable values of a, v, and £ for the magnetic field moving into the 

charge-exchange plasma, R » 1. For this case penetration of the magnetic 

™ 2 23 

field is very slow and the full pressure gradient, B /2p^, can develop, 

Assuming a simple interface, the relation 


2 £_ 
"iV “ 2p, 


( 16 ) 


will give an estimate of where the magnetic field gradient exists. In the frame 

of stationary B, the charge-exchange plasma velocity, v, will roughly be that of 

the solar wind. Solving Eq. (16) for n^, ion density, and using Eq. (1), R is 

found' _o be equal 'vlO^ km, where m. is the charge-exchange ion's mass. Because 

^ -5 

the magnetic field strength is ''■2 x 10 gauss at 1 AU, the Larmor radius is 

3 

'vlO km. The magnetic field will decrease as the radial distance from the sun 

increases and the Larmor radius likewise will Increase. This is a very simple 
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picture for a complex prct)lew. However, it points out that for moat situations 
the apace plasma will exist around the spacecraft and not be affected by the 
solar Mind's magnetic field. 

A-jough measure of the transition distance bet^reen the charge-exchange 
.■jlar.ma and the space i *<as«.»a is where their densa.ties are equal, 'fhc procsidute 
for calculating the Hg charge-exchange plasma has been described £ind has been 
used for calculating the equal-density radius at 90 and 180 deg, to the beam 
direction. This equal-density radius may also be of interest for the ion 
beam. For the latter calculation, the ion beam current was assumed confined 
to a ±15 deg. exhaust cone. The results p.' these calculations are shown in 
Table III. 

The distance 180 deg, to beam direction for Low Orbit is readily available 
in a ground facility. The distance 90 deg. to beam direction^ for low orbit is 
available only at the very largest facilities. All oi her distances eXvjed the 
dimensions of the largest ground facilities. It therefore does not appear 
possible to oimulate this interaction In a ground facility. 

3. Electron Interactions 

As indicated above, c ‘.e electrons are much more mobile than ions, so that 
in a ground vacuuqi facility there will be many electron collisions with the 
facility boundaries for each electron that escapes. In space, there are two 
alternatives. A potential well can be formed by the plasma near tlie spacecraft, 
so that electrons are reflected at the boundaries of this well , similar to 
the reflection at the boundaries of a ground facility. Or the electrons in the 
charge -exchange plasma can escape, and be replaced by ambient space electrons. 


Table III. Equal Density Radii 


Direction Relative 
to Ion Beam 

6 -3 

Low Orbit (10 cm ) 

-3 

Interplanetary {10 cm ) 

0 deg. 

120, w 

40 

km 

90 deg. 

16. m 

5 

km 

180 deg. 

2. m 

0.5 km 
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To show the magnitude of thQ problem, the electron escape rates v^erc calcu- 
lated for tViO case where no reflection occurred. For the Ion beam, tlie. electron 
escape rate was 4i> A. For the entire charge~e>;chnnge plasma It was 34 A. With 
I total electron escape rate of 70 A, and only a 14 A ion escape rate in the 
beam, 4b A of eleo rons v>ouid peed to be either rr '■'lectr I of 1 lected ,1 rom space 
to maintain spacecraft and plasma iieotrality. The rate of electron exchange with 
apace can clearly be quite large. Whether cr not this exchange takes place will 
depend on collision processes for the electrons. 

For interplanetary space, the interaction volume; is quite large (Table III). 

To simplify the calctilations, we will assume a sphere With a radius of 5 km for 

tl»e charge-exchange plasma volume. Integra'* lag th« total number of charge- 

19 

exchange electron/:lon pairs in this volume gives about 1.1 x 10 . At dOOO sec., 

the ion beam adds another 1.5 x lo'^^ electron/lon pairs in this volume. This 

19 

gives a total of about •’.6 v 10 electron/ion pairs within 5 km of the space- 
craft, In comparison, t!ie total miinher of Hg neutrals within this spiitM'e is 
20 

2.1 X 10 . In the 10 - 20 eV electron energy range, the excitaulon/lonlzation 

—20 2 

cro.ss .section of Ug Is several times 10 m , while the Coulomb collision 

-18 2 

cross section of electrons Is of tlie order of 10 m . Clearly, then. Coulomb 
collisions of space electrons with electrons near the spacecraft will be the 
dominant mechanism for capture of t:he.se high velocity, but less dense, space 
electrons. (Coulomb collisions with Ions or elastic collisions with neutrals 
are not significant because of the small energy transferred.) 

At 15 eV, the space electrons have an average velocity of 2.6 x 10^ m/sec. 
7-3 

For a 10 m density, the space electrons should have a collision frequency 
of 

u « nva = 10^ X 2.6 X 30^ x ^ 2.6 x 10“^ (17) 

with each low energy electron in the 3 km radius sphere. For tlie entire 
19 

2.6 X 10 electrons, this would give a capture rate of 

N = 2.6 X 10^^ X 2.6 X 10~^ =» 6.8 x lO^'"' . (IS) 

capt 
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Multiplying this race by the electtonlc charge, q, gives a capture current 
o£ about 100 pA. 

Comparing this electron capture rate with the possible net loss of 56 A 
clearly shows th^’-' capture of ambient space electrons will not be a significant 
proee«?3 in netary space. It seems reasonable to cnnclude > then, t’,.'C 

the two electron populations will -w-oexist in space with negligible i», ;er «ction, 
similar to the two-group theory of electrons in a thruster discharge chamber. 

The low energy electrons must therefore be reflected to maintain charge neutral- 
ity with ions. The ref action won fl be expected near the surface where space- 
craft plasma density equals ambient space plasma denslc-y. 

At low orbit ambient conditions, the situation is considerably different. 

The ambient electron temperature is roughly the same as that of the charge- 
exchange plasma. There should therefore be no significant distinction between 
an escaping charge-exchange electron and a captured ambient space electron. 

The plosma surrounding Che spacecraft should thus blend smoothly with the 
ambient plasma. 

There is another case that is worth considering. A smaller thruster, such 
as might be used for spacecraft attitude control, sUationkeeping, or drag makeup, 
would be expected to have a highar electron temperature. With a higher than 
ambient temperature, the space electrons would be captured while che higher energy 
electrons from the spacecraft would escape. One would therefore expect the 
electron temperature near a spacecraft to be governed by the ambient space 
value when that value is less than the value obtained in a ground test. 

As had been mentioned earlier, the ion (and therefore electron) density is 
related to the Bohm, or ion acoustic velocity. We should address quantitatively 
the question of the effect the ambient plasma may have on this velocity. Two 
plasmas of different density and electron temperature will be occupying the 
same space. What electron temperature will govern the ion velocity? The 
Bohm velocity, 

v/kV^ 
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I’ati be thoiiglit: of aa the vcLof-ity an ion will have after being accelerated 
tlirough a potential i, V^. Therefore, 


v„ = 



( 20 ) 


For the case of a plasma with two Mai"cllinn electron velocity distributions 

liaving different eleet~on temperatures anti densities, tlie following relation 
24 

can bo obtained; 



For typical densities and temperatures of the primary beam plasma and the inter- 
planetary space plasma, the beam ulasma is by far the dominating factor in 
establishing the Bohm velocity. 

For typical values of the solar wind at 1 AU and the charge-exchange plasma 

2 

from a single thruster, the dynamic pressure, ni^n^ v, of the solar wind and 
charge-exchange plasma are roughly equivalent a few tens of meters from 
the spacecraft. This distance will increase as the spacecraft's distance from 
the sun Increases, and us the number of thrusters operating Increases. Because 
one pressure does not overwhelm the other, each individual case of position 
relative to the sun, number of operating thrusters, and solar wind activity will 
have to be calculated to determine this boundary. However, since the mean free 
path for the situations discussed is in the multiple kilometers range, the flow 
of the charge-exchange plasma around the spacecraft should not be affected by 
the solar wind ions. 

4 . Discussion 

The proposed interactions with the space plasma appear reasonable. Ground 
verification, however, ranges from difficult to effectively impossible. Plasma 
calculations have frequently been found invalid in the past. The reason for 
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Cho lack, of validity haa often been collective interactions, which can be very 
diffis'Ult to predict. No collective interactions were considered in this 
analysis. 

F. 'O':;: r, 1, - 

The Hg charge-'‘''change model described appears adequate for preliminary 
estimates of solar array and spacecraft interactions. The electron densities 
predicted fay this model are subject to eieci-roii temperature uncertainties due 
to the variation tbat might be expected with thruster sise, the clustering of 
thrusters and the interaction with the space plasma. There are other effects 
that introduce uncertainty in both electron density and saturation electron 
current density. They are the relative electron and neutral temperatures (which 
result in more or less trajectory bias for the charge exchange ions in the down- 
stream direction) and the departure from axial symmetry for the thruster array. 
All of these effects are estimated to introduce a f.-^ctor of several nneertainty 
in the calculations. In addition, the model is not capable of predicting 
detailed spacecraft shape and surface orientation effects on charge-exchange 
ion trajectories. 

The Hg charge-exchange model, though, is adequate for predicting the 
rough magnitude of tlie electron collection problem at the solar array. It can 
also be used to estimate the plasma environments for spacecraft science instru- 
ments. The arrival rate of Hg ions and possible condensation problems can also 
be evaluated. 

The possibility of surface reactions due to colliding Hg charge-exchange 
ions remains an unknowi. Experience with Hg vapor indicates no problem should be 
expected. Experience v>flth several eV Hg ions, though, is much more Jimited. 

The rate of Hg ion arrival is high enough for much of the spacecraft that cover- 
age could become a serious problem if ions cause reactions on some surfaces 
that Inhibit subsequent vaporization. 

The Mo charge-exchange model has a larger uncertainty than the Hg model. 

It appears, however, that there will be no major problems witii Mo deposition. 

The amount of deposition could still be a problem for certain sensitive instru- 
ments and thermal surfaces near the thrusters. Reducing the uncertainty in Mo 
deposition appears difficult with ground tests. 


To summarize, there appear to be uo major obstacles to the use of ior. 
ilrive. There are areas where spaceeraf t/lon drive liuoractions exi'«t, but can 
be ecntrolleri through the use of proper de.iign approaches. There are also areas 
where uncertainties exist Uiat could cause later problems with sensitive instru- 
ments or surfenc'' PT ./’ot’or of -.‘tee aa^._.tsint.ico in recommendeu <.-at ly in 
the ion drive program to avoid the possible greater expense and schedule impact 
of delayed resolution. 
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